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I .  INTRODUCTION 


T.J.  Karr  put  forward  methods  associated  with  the  use  of  two 
deformable  mirrors  (DM)  to  realize  field  conjugation,  improving 
on  typical  phase  compensation  adaptive  optics  concepts  [ 1 ] . 
Control  algorithms  associated  with  how,  on  the  basis  of  second  DM 
light  strength  distributions,  to  determine  the  phase 
distributions  that  should  be  possessed  by  the  first  DM  are  the 
key  to  this  method.  Karr  presents  “small  signal  algorithms" 
appropriate  for  use  with  uniform  strength  light  sources. 
Basically,  these  only  possess  theoretical  significance.  Actual 
light  sources  very  seldom  possess  uniform  strength.  If  one  wants 
to  make  practical  applications,  it  is  necessary  to  study  more 
universally  applicable  algorithms  associated  with  considerations 
of  light  source  strength  distributions .  What  most  lasers  , 
transmit  are  Gaussian  light  beams.  Therefore,  doing  research  on 
control  algorithms  with  Gaussian  strength  distributions  is  the 
most  practical  thing. 


# 


II.  THEORETICAL  ANALYSIS 


With  a  view  toward  simplicity  and  convenience,  we  discuss 
one  dimensional  problems.  The  methods  are  capable  of  direct 
generalization  to  two  dimensional  forms. 

In  vaccuxims,  near  axis  electromagnetic  wave  motion  equations 


are 


2 


+  v/^=o 


(1) 


Here,  k=2/r/A  is  the  wave  number. 

Let  ,  . 


1  is  wave  length. 


(p=X+ifl> 


(2) 


is  field  logarithmic  amplitudes.  $  is  phase.  Eo  is  a  constant 
electromagnetic  field. 

<p  satisfies: 

2ik  +'V^  <p+(Vx(p)^  =  0 


Acting  as  a  very  rough  estimate,  with  regard  to  standard 
Gaussian  light  beams,  the  expression  for  optical  waist  positions 
<p  is 
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and 


In  order  to  make  approximate  comparisons  with 

^  Tjn  axses,  at  x=y=0  locations,  q>  j  =  6 

At  UI=<Oo/\/T  locations,  are  on  axis  strengths  of 


e-1 


(Vx(p)V(v»  =  i 


with  ir0ga.rd  to  Gaussian  light  hsams/  snsirgies  are  priiriairily 
concentrated  in  the  vicinity  of  light  beam  centers.  Acting  as  a 
preliminary  analysis,  we  first  ignore  ,  making 

linearized  equation: 


lik 


d(p 

dz 


+  Vi*<p  =  0 


(4) 


Taking  <l> 


and  doing  Fourier  transforms: 


I 


+  00 

(p(x)exp(-iKx)dx 

00 


Moreover,  note 


(5) 


Substituting  in  (4),  it  is  possible  to  obtain: 


2ik 


dz 


Solving,  one  gets: 


<p,(z)  =<p^(0)exp(-/ 


K^Z 

~2k 


-K^  (p,=  0 

) 


(6) 


(7) 


Substituting  into  equation  (5),  as  far  as  z  =  z  location 
fields  are  concerned,  it  is  possible  to  ure  z  =  O  location  fields 
expressed  as: 


(piz) 


I 


<;9,(0)exp( 


—  i 


2k 


)exp  (iKx)dK 


(8) 


Comparing  real  parts: 


{[X,(0)cos-|^  +<D^(0)sin-|^]cosKx 


+  [z.(0)sm _  (l>,(0)cos ]  sinKxr 


2k 


(9) 


At  z  =  0  locations,  one  has 


;t(0) 


(0)cosKX-^jO)sinKx}dK 


m^(9). (10)^, 


(10) 


Setting  up  equations  (9)  and  (10)  simultaneously/  it  is 
possible  to  obtain; 


Im[f(x(r))]  -lm[f(x(0))]cos 


Z.(0)=Re[F(z(0))]  + 


2k 


sm 


ZI 

2k 


(11) 


Re(F(z(z))]-Re[F(z(0))]  cos 


O.(0)  =  Im[F(x(0))]  + 


2k 


sm 


2k 


(12) 
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In  the  equations  above/  F  stands  for  Fourier  transform.  Rd  and 
Im  respectively  stand  for  the  real  and  imaginary  parts  of  adopted 
functions.  Because  of  thiS/  from  equation  (5)/  one  knows  that/ 
through  equations  (11)  and  (12)/  phases  on  initial  planes  are 
capable  of  using  exponential  amplitudes  on  transmission  planes 
and  receiving  planes  to  be  represented  as: 

-+« 

O(0)=Im  [z,(0)  +  i‘J^,(0)lexp{iKf}</^  =  Im[F~'(x,(0)+ /<I>,(0))]  (13) 


F-1  stands  for  Fourier  reverse  transforms. 


III.  NUMERICAL  VALUE  SIMULATION  RESULTS  AND  DISCUSSION 


Under  conditions  which  satisfy  r 

formula  (13)  is  universally  applicable.  Here/  numerical  value 
simulations  are  only  carried  out  with  regard  to  Gaussian  light 
beams . 

The  basic  steps  associated  with  numerical  value  simulations 
are  as  follows.  First/  from  the  law  of  conservation  of  energy/ 
one  gets  a  rational  target  reception  field  strength  distribution. 
Second/  from  equation  (13)/  one  constructs  phase  distributions  on 
transmission  planes.  Then/  from  Fresnel  diffractive  integration 
formulae/  calculations  are  done  of  diffractive  strength 
distributions  on  second  DM.  Comparing  differences  between  these 
and  required  fields/  empirical  verification  is  done  of  the 
feasibility  of  algorithms.  Numerical  value  simulation  results 
are  given  below. 

We  adopt  wave  length  A  =  1  micron.  Curve  1  is  transmission 
field  Gaussian  strength  distribution.  Curve  2  is  the  strength 
distribution  required  on  second  DM.  Curve  3  is  a  phase 
distribution  produced  in  association  with  equation  (13).  Curve  4 


is  optical  strengths  obtained  in  association  with  diffraction 
integrals . 

First  of  all,  we  carried  out  numerical  value  simulations 
with  regard  to  required  strengths  having  sine  distribution  forms. 
The  results  were  as  shown  in  Fig.l.  The  produced  strengths  and 
the  required  strengths  are  basically  in  agreement.  Because,  in 
the  case  of  any  curve,  it  is  possible  to  see  forming — in  all 
situations — a  series  of  sine  and  cosine  functions  superimposed  on 
each  other,  this  special  case  possesses  universal  significance. 

In  order  to  avoid  irrational  situations  where  one  has  the 
appearance  of  negative  values  in  association  with  required 
strengths  given,  the  required  strength  distributions  below  are 
given  in  this  way.  From  the  law  of  the  conservation  of  energy, 
one  is  given  a  Gaussian  strength  distribution  with  an  optical 
waist  larger  than  the  optical  waist  of  the  original  light  source. 
Then,  in  the  logarithmic  amplitudes,  sine  and  cosine 
perturbations  are  brought  to  bear. 


x/  m 


Fig.  1  Simulation  on  sine  “shaped  target  intensity  pattern 


n— 0.4m*  0.25m.  r“*  lOOOm 
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Fig. 2  is  nvimerical  value  simulation  results  when  D  =  0.4m, 

<1)0  =  O'. 25m,  and  z  =  1000m.  Due  to  errors  associated  with 
machines  and  algorithms  themselves,  transmission  field  phases 
given  by  equation  (13)  are  made  very  rough.  Moreover,  they  carry 
with  them  small  local  oscillations,  leading  to  abrupt  local 
oscillations'  of  diffraction  strengths. 

What  is  fortunate  is  that  deformable  mirrors  are  very 
smooth.  They  are  capable  of  automatically  getting  rid  of  small 


Fig  .  2  Voilcn.  shake  of  iatensUy  on  .he  sceond  DM  resulted  Irom  .he  slight  local  Huetuafon- 

D*»0.4in.  CD,=  0.25in»  r*  lOOOm 


Fig  .  3  The  intensity  pattern  on  the  second  DM  is  distinctly  improved  after  the  local  phase  ftuctuations 
on  the  first  DM  is  deprived  by  fitting  the  phase  using  spUne  function. 

X)=  0.4m»  O  =  0.25m.  z-  1000m 


(Chinese  page  288  missing) 
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local  phase  perturbations  thereby  improving  diffraction  strength 
distributions.  Deformable  mirror  surfaces  (or  phase 
distribution)  third  (illegible)  degree  spline  functions  are  the 
most  appropriate  for  description.  Making  comparisons  with  Fig. 2 
^19*3  presents  centers  deviated  relatively  more  from  reguired/289 
fields.  Paying  attention  to  increases  in  (Vj.(p)^/V^® 
vslues  along  transverse  distances/  we  estimate  that  tne 
deviations  discussed  above  are  caused  due  to  the  ignoring  of  /r? 
terms  in  equation  (4).  (Vj.<p) 

What  makes  people  pleasantly  surprised  is  that,  although 
dif fraction  strength  distributions  are  very  complicated/  phase 
change  ranges  are  also  relatively  large  (-  8n) .  However/  phase 
curve  forms  are  still  quite  simple.  This  is  extremely 
advantageous  to  deformable  mirrors  accurately  structuring  phases. 
This  is  one  big  advantage  of  Gaussian  light  beams  acting  as  light 
sources. 


Fig. 5  is  numerical  value  simulation  results  when  wo  =  0.15m. 
Comparisons  with  Fig. 3  clearly  show  that  the  larger  D/2wo  iS/  the 
worse  results  are.  The  reason  still  lies  in  ignoring  (VJL<p)^ 
From  comprehensive  consideration  of  the  two  areas  of  energy 
utilization  rates  and  precision/  taking  D/2wo  =  1  is  relatively 
appropriate.  In  this  type  of  situation/  Fig. 6  gives  results  (wo 
=  0 . 2cm) . 
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Fig  .  7  The  shorter  distance  between  the  two  DMs  requires  larger  phase  performance  range  of  the  first  DM. 
Fig.  7 
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Fig. 7  is  simulation  results  when  the  distance  between  the  two  DM 
is  z  =  500m.  Comparisons  with  Fig. 6  clearly  show  that,  taking  z 
and  shortening  it  one  fold,  it  is  then  necessary  to  make 
deformable  mirror  dynamic  ranges  increase  one  fold.  In  actual 
applications,  appropriate  distances  between  the  two  DM  will 
finally  be  determined  by  various  types  of  factors  such  as  target 
strengths  as  well  as  deformable  mirror  operating  ranges,  and  so 
on. 

To  siommarize  what  has  been  said  above,  numerical  value 
simulations  clearly  show  that  this  type  of  algorithm  is 
definitely  capable  of  structuring  appropriate  phase 
distributions,  basically  obtaining  required  target  strength 
distributions.  However,  ignoring  *  introduced  a 

certain  error.  In  particular,  at  places  where  transverse 
distances  are  relatively  large,  errors  are  not  easy  to  ignore. 
Considering  the  influence  of  (Vxq))*  /  it  is  possible  to 

more  accurately  produce  the  needed  strength  distributions.  We 
will  carry  out  further  probes  of  this  later  on. 

Besides  algorithms  needing  further  improvement,  there  exist 
certain  problems  with  plans  themselves.  As  far  as  deformable 
mirrors  are  concerned,  although  capable  of  getting  rid  of  the 
influences  of  small  local  phase  perturbations,  local  strength 
oscillations  created  by  diffraction,  however,  are  still 
unavoidable.  With  regard  to  increasing  deformable  mirror 
dimensions,  it  is  possible  to  make  diffraction  strengths  even 
smoother.  However,  this  will  increase  technical  difficulties  and 
costs.  The  influences  of  local  strength  oscillations  on 


atmospheric  transmission  are  so  large  that  it  is  necessary  to 
make  detailed  probes  of  such  questions  as  whether  or  not  they  can 
be  gotten  rid  of  by  the  wind. 
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